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ABSTRACT 

We present the calibration of the Nuclear Spectroscopic Telescope Array ( NuSTAR ) X-ray satellite. 
We used the Crab as the primary effective area calibrator and constructed a piece-wise linear spline 
function to modify the vignetting response. The achieved residuals for all off-axis angles and ener¬ 
gies, compared to the assumed spectrum, are typically better than ±2% up to 40keV and 5-10% 
above due to limited counting statistics. An empirical adjustment to the theoretical 2D point spread 
function (PSF) was found using several strong point sources, and no increase of the PSF half power 
diameter (HPD) has been observed since the beginning of the mission. We report on the detector gain 
calibration, good to 60 eV for all grades, and discuss the timing capabilities of the observatory, which 
has an absolute timing of ± 3 ms. Finally we present cross-calibration results from two campaigns 
between all the major concurrent X-ray observatories ( Chandra , Swift, Suzaku and XMM-Newton), 
conducted in 2012 and 2013 on the sources 3C 273 and PKS 2155-304, and show that the differences 
in measured flux is within ~10% for all instruments with respect to NuSTAR. 

Subject headings: space vehicles: instruments - X-rays: individual (3C 273) - X-rays: individual 
(PKS 2155-304) 


1. INTRODUCTION 

The Nuclear Spectroscopic Telescope A rray (NuST AR) 
was s uccessfully launched in June 2012 (Harrison et al. 
2013). NuSTA R carries two co-al ig ned co nical Wolter-1 
approximation (Petre & Serlemitos 1985|) Optics Mod¬ 
ules (OMA, and OMB) that focus onto two identical Fo¬ 
cal Plane Modules modules (FPMA and FPMB), each 
composed of four solid state CdZnTe pixel detector ar¬ 
rays (enumerated DET0 through D ET3) wi th a min¬ 
imum detector thres hold of 3keV (|Rana et al. 2009 
Kitaguchi et al.|2011 ). There are 133 shells in each optic. 
The outer 43 shells are coated with a W/Si multilayer 
while the inner 90 shells are coated with Pt/C, limiting 
the highest efficie nt reflective X-ray en ergies below the Pt 
78.4 keV K-edge (Madsen et al. |2009 1. The field of view 
is 13'xl3' and the physical pixel size is 12.3". Sub pixel 
resolution is obtained for events sharing charge among 
multiple pixels, and the physical pixel is subdivided by a 
factor of five in software to an effective pixel size of 2.5". 

Extensive ground calibration and mod eling of thesub- 
systems were performed prior to launch (|Rana et al.|2009[ 


Koglin et al.|20lH|B rcinholt et al.|2011[ Kitaguchi et al. 

2011 Brcjnholt et al. 2012 Westergaard et al. 2012), 


and in this paper we describe the in-orbit refinement of 
the calibration of the NuSTAR observatory. In §[2] we 
present the in-orbit effective area calibration, detailing 
the various components that go into generating the cor¬ 
rect responses; in §[3]we discuss the in-orbit point spread 
function (PSF) calibration; in §[4 we show detector gain 
calibration, and in §[5]we preset the cross-calibration be- 
twee n NuSTAR and the co ncurrent observ atories: Chan¬ 
dra (Weisskopf et al.||2002|), INTEGRAL (Winkler et al. 


et a. 


2003 1, NuST AR ( Harrison et al.|20 13|, Suz aku QMitsuda 


20 07|), Swift (IGehrels et al.| 12004 
NewtorTJ Jansen et ah 2001 H Finally §| 
timing calibration. 


1 ’ ; 
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2. EFFECTIVE AREA CALIBATION 
2.1. The Crab as calibration source 

The Crab is a center filled pulsar wind nebula (PWN) 
powered by a pulsar with a double peaked profile of 
period P ~ 33 ms. It has served as the primary ce¬ 
lestial calibration source for many hard X-ray instru¬ 
ments because of its brightness, relative stability, and 
simp le po wer-la w sp ectrum over the band from 1-100 keV 
(Kirsch et al.|20051. The Crab, however, is too bright for 
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TABLE 1 

NuSTAR Observations log of the Crab 


TABLE 2 
Edges 


Obs ID 

UT Date 

start 

| a Exposure 
(seconds) 

Off-axis angle A/B 
(arcminutes) 

10013022002 

2012-09-20 

2592 

1.52/ 

2.05 

10013022004 

2012-09-21 

2347 

1.51/ 

2.03 

10013022006 

2012-09-21 

2587 

1.49/ 

2.02 

10013023002 

2012-09-21 

2102 

4.01/ 

3.38 

10013024002 

2012-09-20 

2258 

4.77/ 

4.20 

10013025002 

2012-09-25 

1235 

5.08/ 

5.68 

10013025004 

2012-09-25 

1161 

5.02/ 

5.62 

10013025006 

2012-09-27 

1593 

4.95/ 

5.55 

10013026002 

2012-09-26 

2540 

5.56/ 

5.64 

10013026004 

2012-09-26 

1162 

5.57/ 

5.65 

10013027002 

2012-09-27 

1182 

5.43/ 

4.99 

10013027004 

2012-09-27 

1105 

5.48/ 

5.03 

10013028002 

2012-09-28 

1601 

6.30/ 

5.82 

10013028004 

2012-09-28 

1254 

6.31/ 

5.82 

10013029001 

2012-10-25 

2909 

4.07/ 

3.41 

10013030001 

2012-10-25 

3023 

4.65/ 

5.29 

10013031002 

2012-10-25 

2507 

2.19/ 

2.86 

10013032002 

2012-11-04 

2595 

1.09/ 

1.27 

10013033002 

2012-12-19 

1383 

1.30/ 

1.93 

10013033004 

2012-12-21 

1269 

1.19/ 

0.61 

10013034002 

2013-02-14 

988 

1.02/ 

1.69 

10013034004 

2013-02-14 

5720 

0.54/ 

0.88 

10013034005 

2013-02-15 

5968 

0.75/ 

0.87 

10013035002 

2013-02-15 

9401 

5.32/ 

5.44 

10013037002 

2013-04-03 

2679 

2.14/ 

2.65 

10013037004 

2013-04-04 

2796 

3.28/ 

3.95 

10013037006 

2013-04-05 

2799 

3.24/ 

3.90 

10013037008 

2013-04-18 

2814 

3.28/ 

3.95 

10013038002 

2013-04-08 

3084 

4.20/ 

4.78 

10013038004 

2013-04-09 

2217 

4.21/ 

4.79 

10013038006 

2013-04-10 

266 

4.42/ 

5.01 

10013038008 

2013-04-17 

2231 

4.16/ 

4.74 

10013039001 

2013-05-01 

5206 

3.79/ 

3.39 

10013039002 

2013-05-02 

590 

3.66/ 

3.25 

10013039003 

2013-05-03 

583 

3.66/ 

3.25 

80001022002 

2013-03-09 

3917 

1.42/ 

1.75 

10002001002 

2013-09-02 

2608 

1.72/ 

2.09 

10002001004 

2013-09-03 

2386 

1.73/ 

2.26 

10002001006 

2013-11-11 

14260 

1.08/ 

1.28 


a Effective exposure time corrected for dead-time. 


reliable pileup corrections for most CCD based focusing 
X-ray instruments, and has been replaced with fainter 


sources, such as the Crab-like PWN G21.5 +0.9 ( Tsu.ji- 
mo to et al.|2011 ), the blazar PKS 2155-304 (Ishida et al. 


201ip or the quasar 3C 273. The latter two are variable 


and only serve for cross-calibration among observatories, 
while G21.5+0.9 would require unfeasible long integra¬ 
tion times for NuSTAR to obtain the required statistics 
to calibrate the response at 80keV. The Crab therefore 
remains the best choice for the internal calibration in the 
X-ray band covered by NuSTAR. 

The spatially phase-averaged integrated spectrum of 
the Crab nebula+pulsar in the ITOOkeV X-ray band 
has been well-described by a power-law with photon 
index of T ~ 2.1 (RXTE , BeypoSAX , EX OS AT, IN- 
TEGRAL/ JEM_X; |Kirsch et al.[|2005|). Above 100keV 
the hard X-ray instruments [IN 1 E (JR .4 L/SPI/1S G III, 
CGRO) measure a softer index of T ~2.20-2.25, and be¬ 
low lOkeV the instruments with CCD detectors a harder 
spectrum (P < 2.1). The softening above 100keV pri¬ 
marily comes from the curvature of the pulsed spec¬ 


Edge (keV) 

W 

Pt 

^3 

10.21 

11.56 

Li 

11.54 

13.27 

Li 

12.10 

13.88 

K 

69.53 

78.40 


trum (Kuiper et al. 2001). The hardening in CCD X- 
ray instruments may be due to photon pile-up and al¬ 
though models exist for these instruments to deal with 
the pile-up the Crab usually still challenges the pile-up 
cor rections and re quires special non-standard reductions. 
Weisskopf et al. (2010) did a study with RXTE/ PCA, 
XMM, and Ai'OAlooking for deviations from a power- 
law in the 0.2-50 keV energy range and concluded that 
within the precision of the available instrumentation 
there is no detectable bend in the phase-averaged inte¬ 
grated spectrum of the Crab. Over the 16 years RXTE 
was operational and regularly monitoring the Crab, the 
spectral index was see n to vary by a p eak- to-p eak vari¬ 
ation of Ar ~ 0.025 (Shaposhnikov et al. 2012). This 
variation is consistent with the observed spread between 
instruments, but it is slow and on average over the 16 
years the Crab has remained at T = 2.1. 

Because the non piled-up instruments covering the 1— 
100 keV band agree on a photon index of T = 2.1 ± 0.02 
and do not measure any curvature in the spectrum across 
this band, we have calibrated the effective area against 
a Crab index of T = 2.1. We model the differential pho¬ 
ton spectral model of the phase-averaged integrated spec¬ 
trum of pulsar+nebula as a simple absorbed power law, 


dN(E) 

dE 


= tbabs (E)NE~ l phs _1 cm _2 keV~ 


( 1 ) 


Here E is the photon energy, tbabs is the inters tellar 
absorption using Wilms abund ances j Wilms et al.pOOO) 
and Verner cross-sections (Verner et al. 1996j), T is the 
power-law photon index, and N is the normalization fac¬ 
tor. When r = 2.1 and N = 10 keV -1 cm -2 s -1 at 1 
keV, we refer to this as the canonical Crab model. 

In terms of flux stability the Crab has been tracked 
in great detail over the last decade. It remain s largely 
stable, albeit with flux changes as much as 7% (Wilson- 
Hodge et al.)2011) across the 10-100 keV bandpass. This 
corresponds to a decay in the flux of 3.5% per year over 
the period it has been observed. The NuSTAR observa¬ 
tions span approximately a year, bu t du e to the optical 
axis knowledge (to be discussed in §2.4) we expect flux 


differences between repeated observations to be on the 
same order (5%), and therefore the intrinsic variability 
of the Crab is not an issue. In addition, after adjust¬ 
ment of the effective areas of the two mirror modules as¬ 
suming the canonical Crab model, we adjust the overall 
normalization by 15% and use the two targets PKS 2155- 
304 and 3C 273 observed simultaneously with Chandra, 
NuSTAR , Swift ., Suzaku and XMM-Newton to confirm 
the final global normalization of our responses with re¬ 
spect to the other observatories. 


2.2. Observations and Data Reduction 

The NuSTAR Crab calibration campaign measured 
the Crab spectrum at different locations on the detec- 
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Fig. 1. — Ratio of the base effective area to model (r=2.1, N=10 
ph keV -1 cm -2 s -1 ) for observation 10002001006, which was the 
longest single continuous observation taken at a fixed off-axis angle. 


tor modules at various off-axis angles. This resulted in 
the 39 individual observations listed in Table [I] 

Despite the angular extent of the Crab of ~120" x 100", 
we treat it as a point source since the center of the neb¬ 
ula where the pulsar resides dominates the emission. The 
difference in the effective area response between using a 
point source effective area centered on the pulsar and a 
count-rate adjusted extended effective area is only ~1- 
2%. This difference is largely constant as a function of 
energy and is therefore an acceptable error. Treating 
the Crab as a point source allows us to use the built- 
in pipeline corrections of the Aperture Stop (AS), Ghost 
Ray (GR) and PSF correction s (the details of these com¬ 
ponents will be discussed in §2.3.2). These corrections 
are not applicable to extended source responses, and 
since they are much larger than the error introduced by 
assuming the Crab is a point source, they are extremely 
important to include. 

We extracted counts from the largest possible source 
region, typically 200", that did not extend beyond the 
edge of the detector. At this radius about 95% of all pho¬ 
tons are contained in the source extraction region. We 
used the NuSTARDAS pipeline version vl.3.0 to reduce 
the data, generate the detector response, RMF, and to 
apply GR, AS and PSF corrections to the base effective 
area. 


2.3. Analysis 

Prior to launch, we used ground calibration and mod¬ 
eling to produce a physic s-based optics ray-trace model 
( Westergaard et al.[2012[ ) to generate a set of base effec¬ 
tive area and PSF files. Figure [l] shows the ratio of the 
canonical Crab model folded through the base responses 
to one of the observations. This comparison shows sev¬ 
eral features that indicate inaccuracies in the ground cal¬ 
ibration: 


10 keV and 14 keV the Pt and W L-edge complex, 
along with the top few coatings of the multilayer 
stack, create residuals. 

• Below 20keV X-rays reflect off the optics by total 
external reflection. Above this energy the multi¬ 
layer coatings become effective and the resonant 
features, enhancing the reflectivity, are sensitive to 
the structure of the multilayer that we can only 
approximate because the coatings are not precisely 
tabulated for each individual mirror segment of 
which there are more than 6000. 

• Between 30-40 keV the effective area deviates grad¬ 
ually from the model areas. The exact reason for 
this discrepancy is not known. 

• At 69.4 keV the W K-edge creates a residual for 
the same reasons as the L-edge. At 78.4 keV the 
Pt K-edge is just evident. 

In the following sections we model the detected counts 
in a given instrumental pulse height bin, C(PI,9), ac¬ 
cording to the equation, 

C(PI , 9)=J ^j^R(PI, E, 6)dE. (2) 

Here dN(E)/dE is the model differential photon spec¬ 
trum of the observed target as a function of incident 
photon energy, E (flux in units of photons s -1 cm~ 2 ), 
and R(PI, E, 6) is the response matrix that captures the 
photon indicent in a given pulse height bin, PI (units of 
cm 2 ), at an off-axis angle, 9. In practice, this integral 
is approximated as a finite sum by sampling R(PI, E, 9) 
on a fine grid. The off-axis angle, 9 , is the angle of in¬ 
coming X-rays with respect to the optical axis of the 
mirrors. As the optical axis moves with respect to the 
detector position during an observation, the modeled re¬ 
sponse is sampled on a finite time grid and then summed 
for a given exposure. 

As is typical for X-ray astronomy missions, the re¬ 
sponse matrix is divided into two components, 

R(PI, E, 9) = RMF(PJ, E)A{E, 9). (3) 

RMF(P7, E) is known as the redistribution matrix, 
which contains detector quantum efficiency and resolu¬ 
tion effects, and is unitless (a fraction between 0 and 1, 
where a value of 1 indicates 100% quantum efficiency). 
The quantity A(E, 9) is the effective area, also known as 
the ancillary response function (ARF), which captures 
the effective area of the mirror optics, as well as several 
other factors unique to NuSTAR, 

A(E , 9) = A 0 (E, 9)V(E, 6»)detabs(P)GR(P, 9) 

AS (E,0)C(E,0). (4) 


• Below 10 keV the response turns over exponentially 
due to uncertainties in the thickness of the Plat¬ 
inum (Pt) electrode on the detector surface. 

• At around 10 keV the Tungsten (W) L-edge causes 
residuals (See Table [2). The exact shape of the 
W edge depends on the optical constants of the 
multilayer coatings, which are very sensitive to the 
crystalline structure and density of the W. Between 


Here Aq is the modeled effective area of the mirror seg¬ 
ments estimated from the physics-based ray-tracing sim¬ 
ulations, and V is the geometric vignetting function, also 
based on the ray-tracing simulations. The quantities 
detabs, GR, and AS are the detector dead layer absorp¬ 
tion, ghost ray correction, and aperture stop correction, 
respectively, described in the next sections. C is an em¬ 
pirically deri ved co rrection factor, which we will discuss 
at length in § |2 3.3 
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2.3.1. Detector Absorption Correction 

A Pt contact coating and a CdZnTe dead layer on the 
detector surface cause absorption at low energies. No 
ground calibration of the absolute efficiency at these en¬ 
ergies was performed due to cost and scheduling con¬ 
straints. Therefore uncertainties in the low energy de¬ 
tector efficiency remained prior to launch. To investigate 
this detector absorption, we built an XSPEC absorp¬ 
tion model with cross-sections created by Geant4. The 
adopted photon interaction model is the Livermore low- 


library, EPDL97 ( 

Cirrone et al. 

2010). 

Using the Crab complicates tl 

re calibration of this ab 


sorption layer since the Ah column creates a similar ef¬ 
fect. At an JVh ~ 3 x 10 21 cm -2 the contribution between 
3-4 keV is significant enough that to derive the detector 
absorption ( detabs ) parameters we had to use 3C273, 
which has a low Ah of 1.79 x 10 20 cm -2 (Dickey & Lock- 


man 1990j) . At these Ah values the quasar 3C 273 has no 
measurable absorption above 3 keV. We fitted only from 
3 to 9 keV, avoiding the W L-edge, and since reflection off 
the optics at these energies occurs in the regime of total 
external reflection, where the exact composition of the 
multilayer stack does not matter, the ground calibration 
responses were acc urate en ough for our purpose. 


We used XSPEC (Arnaud 1996) for the analysis, fitting 
an absorbed powerdaw to the data, employing Wilms 


abundances (Wilms et al. 20001 and Verner cross-sections 


(Verner et al. 1996J). Using a power-law index of T = 
: .63, derived from the cross-calibration campaign to be 
described in 1 5j we fitted the thickness’s of the Pt and 
CdZnTe for DETO as shown in table [3] for FPMA and 
FPMB. Because we only have observations with 3C273 
on DETO, we used these numbers for all the Crab obser¬ 
vations made on DETO to derive an Ah = 2.224 x 10 21 
cm -2 for the Crab. We applied this Ah to simultaneous 
fits of all the Crab observations on each detector respec¬ 
tively for each FPM to derive the other detector depen¬ 
dent absorption parameters. These are summarized in 
Tabic H 

The measured absorption coefficients are formed into 
energy-dependent correction files, which are multiplied 
onto the effective area as an absorption layer in the opti¬ 
cal path. These energy correction files, called detabs, are 
stored in the NuSTAR calibration database (CALDB) 
maintained at NASA HEASARCQ 


2.3.2. Ghost Rays and Aperture Correction 


TABLE 3 
detabs PARAMETERS 


Module 

Detector 

Pt ((im) 

CZT H 

A 

0 

0.1143 

0.2275 

A 

1 

0.0993 

0.3069 

A 

2 

0.0712 

0.3476 

A 

3 

0.1081 

0.2554 

B 

0 

0.1209 

0.2451 

B 

1 

0.0867 

0.2887 

B 

2 

0.1091 

0.2430 

B 

3 

0.0869 

0.3248 


1 https://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/nustar/ 
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Fig. 2.— Aperture stop correction, AS (E, 9). The correction is a 
function of azimuthal location of the source with respect to the op¬ 
tics and off-axis angle. The spectral dependence occurs because the 
aperture stop preferentially blocks the low energy photons emerg¬ 
ing from the outer shells of the optic, which have a larger area. 
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The NuSTAR observatory has two benches, the Focal 
Plane Bench (FB) and the Optics Bench (OB) separated 
by 10.14 meters, which move relative to each other due 
to motions of the mast that connects them and causes 
the optical axis to travel across the FPMs. Because the 
optics act as thin lenses, this track is not equivalent to 
the path made by the source; the rotational component 
of a bend in the mast does not project into a movement 
of the source on the detectors, only the translational part 
does. In the system there are therefore two paths that 
must be kept track of as a function of time: the optical 
axis (OA) path, which is typically on the order of 30"- 
60" wide, and the source path, which can be considerably 
larger and is mainly driven by the pointing stability of 
the spacecraft bus. The calculation of off-axis angle is 
carried out as a function of time, and since the resolution 
of the response files is 10", the off-axis aspect histogram 
is binned at this resolution. 

Three aperture stops are attached to each FPM with 
the top acting as the limiting aperture. The diameter of 
the aperture opening is 58 mm and the total thickness of 
the solid part of the aperture is 2 mm, layered into 0.75 
mm Al, 0.15 mm Cu and 1.1 mm Sn. The purpose is to 
block the diffuse background by limiting the field of view 
(FOV) of the detectors to the open sky. However, this 
comes at the cost of a reduction in effective area at off- 
axis angles due to the clipping of the edges of the optical 
path between detector and optic. Since the aperture stop 
is fixed with respect to the detector, and the OA is mov¬ 
ing independently, the correction is strongly observation 
dependent. This clipping is a purely geometrical effect 
and can be calculated by a ray-trace. We have assembled 
the result as a function of off-axis angle and distance to 
aperture stop center in a look up table used by the soft¬ 
ware that creates the ancillary response file for a given 
observation (numkarf). An example of the magnitude of 
this correction is shown in Figure [2] The correction is 
azimuthally dependent since the aperture stop and the 
OA are not centered, and the correction has a spectral de¬ 
pendence, preferentially cutting away low energy photons 
since the majority of these come from the outer shells of 
the optic and thus are more prone to being blocked than 
light focused by the inner shells. Figure [2] shows that the 
correction is only important for off-axis angles > 3 7 . 
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Fig. 3. — Ray-trace Ghost Ray (GR) patterns. Images shows the distribution of GR as a function of off-axis angle. The source location 
on the detector is initially on-axis and hits the center of the detector marked by the black cross. As the source moves off-axis towards 
the left the GR population moves towards the right. The peak intensity of the GR initially moves towards the right, but as soon as the 
source and the GR population meet, the peak intensity of the GR moves together with the source out of the detector’s field of view. The 
petal-like pattern continues to grow and becomes more diffuse as the source moves further off-axis, until at around 40' reflection of GR is 
no longer possible. 


The NuSTAR optics are based on a Wolter-I conical 
approximation, which is a grazing incidence, double mir¬ 
ror system. A properly focused photon will reflect twice 
off the optics before exiting. However, it is possible for 
a photon to reflect only once by either the upper conical 
section or the lower conical section. This occurs at either 
very shallow or very steep angle and we call these single¬ 
bounce photons Ghost Rays (GR). The pattern of GR is 
very distinct, as shown in Figure [3j and axi-symmetric. 
Figure [3] shows the pattern for a series of off-axis angles, 
starting with the source being on-axis and hitting the 
center of the detector. In this case the GR are uniformly 
distributed in a circular halo, which we never see in NuS¬ 
TAR because it is blocked by the aperture stop. As the 
source moves off-axis and the spot wanders towards the 
left side of the detector, the GR can be observed entering 
from the left until around 6' they overlap with the double 
bounce photons. Since the GR are now inside the typical 
extraction region, the net effect is an increased effective 
area. Like the properly reflected photons, these are also 
subject to being partially blocked by the aperture stop, 
and so this must also be taken into account. Due to the 
complexity of the GR pattern, the corrections were only 
derived for circular regions. Figure [4] shows an example 
of the correction, which, like the AS correction, is a func¬ 
tion of azimuth, but also depends on source extraction 
region size and energy. The suppression of the low en¬ 
ergies is because of the aperture stop correction to the 



Fig. 4.— GR correction at 4'. The correction is a function of az¬ 
imuthal angle of the source with respect to the optics, off-axis angle 
and extraction region size. It only works with circular extraction 
regions. The correction increases the effective area. 

GR. 

The energy dependent correction is applied to the an¬ 
cillary file in numkarf, and because of the complexity of 
the GR correction files and their dependency on extrac¬ 
tion region, this correction is only available for circular 
and annular regions, and only for point sources. 

2.3.3. Empirical Correction and Vignetting 

After adjusting for all the above terms ( detabs , AS, 
GR), there still remain systematic residuals in the Crab 
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Fig. 5. — Effective area correction factor, C(E,9), as a function 
of off-axis angle and energy for FPMA (top) and FPMB (bottom). 
Values are shown for O' (solid), 3' (dotted) and T (dashed). 

spectrum which are not accounted for, with typical resid¬ 
uals of a few percent for on axis targets and 10-20% for 
off axis targets. 

Although it is desirable to have a completely physics- 
based response matrix derived from a priori information, 
this is not usually achievable with limited project re¬ 
sources. We found that the residual trends as a function 
of energy, E, and off axis angle, 9 , are very reproducible, 
and therefore we elected to adjust the response by an 
empirical correction factor C ( E , 9 ), using the canonical 
Crab spectral model. C is a fractional quantity, where 
a value of 1 indicates no correction. We note that C is 
a function of energy E (not pulse height), so its effect 
is equivalent to modifying the effective area term in the 
ancillary response function (ARF), and its effects must 
be folded through the response matrix to be seen in the 
model pulse height spectrum. As a purely practical mat¬ 
ter, C(E , 9) is stored along with the vignetting function 
in the NuSTAR CALDB (i.e., a new vignetting function, 
V'(E,9) = V(E,9)C(E,9)). 

We investigated several functions for C(E, 9) and de¬ 
cided to use a piece-wise linear interpolation function, in 
both energy and off-axis angle. This function consists of 
piece-wise linear segments between fixed control points 
E C i in energy. We chose the control points of the spline 
based on the known features. Figure [5] marks the loca¬ 
tions of the nodes with 27 control points between 3.0 and 
80 keV. The function is completely specified by the value 
at each of the control points, which are considered the 
parameters of the correction model. The angular depen¬ 
dence is treated as linear as well, as described further be¬ 
low. We investigated using Chebyshev polynomials or cu¬ 
bic splines as an interpolation function, but encountered 
difficulty with ringing and other unconstrained behavior 
in the neighborhood of the control points. The interpo¬ 
lation approximation results shown here worked well for 
our purposes to within available statistics of the data. 
An improved model may be found in the future, but the 




Fig. 6.— Effective area correction factor, C(E,6 ), derived from 
individual Crab spectra, at representative energies: 3 keV (top), 
10.1 keV (middle) and 46.3 keV (bottom). FPMA is black and 
FPMB is red. The thick lines represent the best fit interpolation. 
Error bars are standard errors from individual Crab fits, and ap¬ 
pear to underestimate the true scatter. Some variability, such as 
the few percent dip around 3 - -4 arcmin at lower energies, is not 
captured by the simple linear model. 

present correction achieves residuals at the few percent 
level over almost the entire NuSTAR energy range. 

We derived the correction spectrum in two stages. 
First, we fit individual Crab spectra to derive the energy 
dependence for several off-axis angles. Second, we used 
linear interpolation to model the angular dependence to 
derive a single global correction model. 

For the first stage of analysis, we combined the ob¬ 
servations in V bins to improve statistics, resulting in 
seven spectra between 0' (on-axis) and 7' off-axis. We fit 
the parameters of the correction model while holding the 
canonical spectral model fixed, using standard \ 2 fitting, 
and requiring at least 200 counts per bin. The result is a 
linear model for each of the seven spectra. We note that 
for the highest control points at 77 and 80 keV, there 
were not enough counts to reliably estimate the param¬ 
eter values at both points independently, so the value of 
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Energy (keV) 

Fig. 7. — Effective Area curves (top) for OMA, and OMB, 
A(E,9), for CALDB file version v006. Vignetting curves (bot¬ 
tom), V'(E,0), of OMA and OMB for off-axis angles 1', 2' and W 
for CALDB file version v006. 



Fig. 8.— Ratio of the base effective area to model (r=2.1, N=10 
ph keV~ 1 crrr 2 s" 1 ) for all Crab observations within 3b 

the response at the 80 keV control point was fixed manu¬ 
ally to be 90% of the value at 77 keV. In the final result, 
we claim calibrated flux levels only between 3-78 keV. 

The corrections show an angular dependence nearly 
linear with off-axis angle (Figure [6]). This indicates that 
the physics-based ray-tracing captured most of the angu¬ 
lar dependence of the reflectivity of the mirrors, except 
for a residual trend. We therefore treat the angular de¬ 
pendence of the correction as linear between on-axis and 
7' off-axis. For off-axis angles above 7' (the outer parts 
of the focal plane) we lack the Crab data to constrain 
the model, so we cannot verify if the linear trend contin¬ 
ues. To be conservative, we do not extrapolate the linear 
trend beyond 7', but simply hold the function fixed at 
the 7' value. 


TABLE 4 

Measurement Repeatability Statistics 


Subset 

Mean 

Std. Dev. 

Relative Flux 

All 

-0.4% 

2.6% 

On-Axis (6 <3’) 

-0.9% 

3.0% 

Off-Axis (0 >3’) 

-0.9% 

3.0% 

DET0 

-0.3% 

2.5% 

DET1 

-0.7% 

2.5% 

DET2 

+0.8% 

1.4% 

DET3 

-0.0% 

3.4% 

DETO+3 Straddle 

-1.7% 

4.0% 

FPMA-FPMB Diff. 

-0.0% 

2.6% 

Power-Law Photon Index, F 

All 

— 

0.009 

FPMA-FPMB Diff. 

— 

0.012 

Absorption, A r H(10 22 cm 

-*) 

All 

- 

0.150 


In a second stage of analysis, we fitted the individ¬ 
ual correction values at each control point independently 
with a linear function, using the standard errors from 
the individual fits, to derive one global correction func¬ 
tion for each FPM. We stress that this global correction 
is an ensemble-derived model using all of the Crab obser¬ 
vations; no single Crab spectrum was used to fix the final 
global model parameters. Thus, the final global correc¬ 
tion model consists of 28 correction parameters on-axis 
and 28 parameters 7' off-axis, for a total number of model 
parameters of 56 for each FPM. These are derived from 
7 x 26 (=182) individual correction spectra, which were 
in turn derived from 7 x 940 (~ 6600) individual Crab 
count bins for each FPM. Figure 15] shows the correction 
function C(E,9) as a function ofenergy E and off-axis 
angle 9 for both FPMs. The systematic nature of the 
adjustments and some of the features mentioned in § |2.3| 
are evident. 

As noted above, the correction factor C(E,9) and 
geometric vignetting function derived from ray tracing 
V(E,9) were combined into a single effective vignetting 
function V'(E,9), and stored in the NuSTAR CALDB. 
This effective vignetting function is shown in Figure [7] 
(bottom), along with the resulting total effective area 
for each module (top). 

We reprocessed the Crab data with the adjusted cali¬ 
bration files to verify that the correction was successful. 
Figure [Si shows the residual spectrum for all Crab data 
within T of on-axis for both modules. The residuals are 
typically better than ±2% up to ~40 keV. Between 40- 
80 keV, residuals are dominated by counting statistics of 
the Crab data, but are typically 5-10%, so the systematic 
errors are less than that. 

The final correction achieves few percent residuals over 
the entire NuSTAR energy range of 3-78 keV, and is 
applicable for off-axis angles of 0' 7'. We do not advise 
to fit spectra above 78 keV. 

2.4. Measurement Repeatability 

Having calculated the correction function, as a consis¬ 
tency check we reprocessed all individual Crab data sets 
to determine if the specified Crab spectral parameters 
were recovered. Although the Crab spectral shape and 
normalization were enforced at the ensemble level during 
the correction process, we did not required that individ- 
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Fig. 9. — Residual spectral parameter values for relative flux 
(top) and power law photon index (bottom) as a function of off- 
axis angle, 0. FPMA is black and FPMB is red. The canonical 
values are shown as dashed lines. 


was stable to about ±0.01 (lex); this can be taken as an 
approximate estimate of repeatibility errors in spectral 
slope. For neutral absorption, NuSTAR is not particu¬ 
larly sensitive to the Crab’s hydrogen column; the sample 
standard deviation is comparable to the canonical value 
itself. 

The fluctuations in flux and slope are understood, and 
are consistent with the estimated error in optical axis 
location. The knowledge of the optical axis location is 
accurate to ~ 30", and because the optical axis is con¬ 
stantly moving due to the mast motions, there will be 
situations where the agreement between effective areas 
are better than in others. The flux level is also affected 
by the understanding of the PSF (to be discussed in the 
next section), and the fidelity of the instrument map; 
gaps, dead pixels etc. The gap width between detectors 
was measured on the ground by scanning an X-ray beam 
across the gaps, but due to systematic errors getting the 
exact physical gap width is challenging, making it dif¬ 
ficult to correct the flux exactly in the cases where the 
core of the PSF falls, or crosses the gap. 

Observed differences between the two modules of up 
to 5% are not uncommon and are simply a result of the 
optical axis configuration during the time of observation. 
In extension, observed flux differences less than 5% be¬ 
tween different observations of the same target, should 
not be interpreted as intrinsic to the source. 


ual spectra conform, nor did we require that FPMA and 
FPMB produce the same results. Any deviations at the 
individual observation level would indicate unmodeled 
systematic errors, as well as errors in the assumption of 
a constant Crab. To estimate this, we computed sample 
biases and standard deviations for various subsets of ob¬ 
servations. This is a measure of how repeatable NuSTAR 
observations are. 

We fit each of the spectra with the canonical spectral 
shape, but allowed the parameters to vary. The param¬ 
eters were the relative power-law flux normalization and 
the power law photon index, T, while the neutral absorp¬ 
tion parameter, Ah, was fixed. We also performed a sep¬ 
arate fit where in addition Ah was allowed to vary. Re¬ 
sults from FPMA and FPMB were recorded separately. 

We recorded the individual parameters and computed 
statistics of various subsets of observations. The results 
are shown in Table [4j In addition to the sample varia¬ 
tion among “All” observations, we also examined subsets 
where the target was near on-axis and far off-axis, and on 
individual detectors (DETO-3). These designations are 
not mutually exclusive. In a few cases the Crab strad¬ 
dled two detectors, and these are shown separately. We 
also computed the sample variations of the differences 
between FPMA and FPMB. For relative flux, we com¬ 
puted the sample bias as well as the standard deviation; 
for photon index and absorption, we only list the sample 
standard deviation. The results are also shown graphi¬ 
cally in Figure 9] as a function of off-axis angle. 

The results demonstrate that average relative flux bi¬ 
ases are small, typically <1%. Scatter in the relative 
fluxes is typically 2-3% (1 <t). There are no large flux 
offsets between FPMA and FPMB during the same ob¬ 
servation; in the worst case the difference was 3.8% in 
this sample. There are no strong biases as a function of 
off-axis angle or detector number. The power law index 


3. POINT SPREAD FUNCTION CALIBRATION 
3.1. Near On-Axis PSF Shape 

Calculating the correct effective area requires a good 
understanding of the PSF. Prior to launch a ray-trace 


TABLE 5 

NuSTARPSF Observations log 


Target 

Obs ID 

Exposure 

(ks) 

Off-axis anglaq 
(') 

Cyg X-l 

10002003001 

9 

0.5 

GRS 1915+105 

10002004001 

15 

1 

Vela X-l 

10002007001 

11 

1 

Vela X-l 

30002007002 

7 

1 

Vela X-l 

30002007003 

24 

2 

GS0834-430 

10002018001 

31 

1 

Her X-l 

30002006002 

28 

2 

Her X-l 

30002006005 

22 

1 

Her X-l 

30002006007 

27 

1 

Cyg X- f\ 

00001007001 

2.4 

3 

Cyg X-1 V 

00001008001 

4.3 

3 


a Average value for the distribution 
b for off-axis PSF 


TABLE 6 

Measured half power diameter for the on-axis 
OBSERVATIONS IN TABLE [5] 


Energy 


FPMA 


FPMB 

(keV) 

HDP 

(") 

CTHPD 

(") 

min. 

n 

max. 

(") 

HPD 

(") 

O'HPD 

(") 

min. 

(") 

max. 

(") 

3-4.5 

70.3 

2.4 

66.7 

75.5 

65.6 

2.4 

62.9 

69.6 

4.5-6 

67.1 

1.0 

64.7 

67.7 

62.6 

1.2 

60.9 

64.7 

6-8 

64.7 

1.0 

62.8 

65.7 

60.7 

1.4 

58.8 

63.7 

8-12 

63.5 

1.1 

61.8 

64.7 

59.5 

1.5 

57.9 

62.8 

12-20 

63.4 

1.1 

61.8 

64.7 

60.3 

1.2 

58.8 

62.8 

20-79 

63.4 

1.0 

61.8 

64.7 

62.4 

1.5 

60.8 

65.7 
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model of the PSF was developed and tested against lim¬ 
ited grou nd calibration scatter data from individual flight 
mirrors (Westergaard et al.| 2012). This original set of 
calibration hies formed the basis of pre-launch PSFs de¬ 
termined every 0.5' out to an off-axis angle of 9 = 8.5'. 

To conduct the in-orbit calibration of the PSF, we ob¬ 
served the strong point sources listed in Table [5j The 
main point of these observations was to establish the cor¬ 
rections required to the ground based estimates. To this 
end, we extracted events from the point source observa¬ 
tions in several energy ranges (3-4.5, 4.5-6, 6-8, 8-12, 
12-20, and 20-78 keV) and produced radial profiles cen¬ 
tered at the peak of the intensity map. Table [6] shows 
the half power diameters (HPD) for the nine on-axis ob¬ 
servations. 

We compared these radial profiles with the ray-trace 
model. Since the off-axis angle changes during each 
observation due to the mast motions, we weighted the 
ray-traced PSF with the aspect solution before compar¬ 
ing. We fitted the observed radial profiles to the aspect- 
weighted ray-trace PSF and found that: (i) the central 
cores of observed radial profiles are slightly broader than 
that of the model PSF, (ii) there is an unaccounted for 
wing at R > 100" in the radial profiles. The broaden¬ 
ing of the central core does not affect the FWHM and is 
an effect of incomplete astrometric aspect reconstruction. 
The origin of the wider wing is unknown. 

To take care of the central broadening we convolved 
the modeled PSF with a Gaussian function. Since the 
origin of the wing is unknown, we decided to empirically 
model it by introducing a component on top of the ray- 
trace PSF. Because background subtraction is not feasi¬ 
ble during the fitting, we had to include a background 
component as well, which is well known and can be mod¬ 
eled accurately ( Wik et ahj2 014). The backgrounds were 
in general smali compared to the sources, dominating 
only above R > 500". 

The 2-D model PSF can be expressed by the following 
formula: 


y) = C'o'P RT (z, y) ® G(x, y, a) 
+C ie - a ^ x2+y2 + C 2 B{x, y) 

where Trt is the ray-trace PSF, G(x,y,a) is the Gaus¬ 
sian function, e a ^ x2+y2 is the wing component, B{x,y) 
is the background model, and Co, 1,2 are normalization 
constants. 

To fit the radial PSF profiles, we first convolved the 
ray-trace PSFs with the Gaussian function and weighted 
them using the aspect solution of each observation to 
produce an aspect-weighted PSF. The wing model was 
added to construct an effective PSF, and we multiplied 
it with the exposure map for the observation. Since 
the point source exposure map is almost flat, the over¬ 
all shape of the PSF does not change much. However, 
multiplying with the exposure map is important for some 
localized effects such as detector gaps. 

The fitting parameters are the amplitude of the ray- 
trace PSF (C 0 ), width of the Gaussian convolution core 
(cr), amplitude and decay constant for the exponential 
wing (Ci and a), and amplitude of the background (C 2 ). 
Since fitting the observations simultaneously was diffi¬ 
cult given different aspect histories (i.e., different effec¬ 



Fig. 10.— Radial profile of CygX-1 in the 5-8keV band. Fit 
model components are also shown. Background components were 
fit separately but only the combined model is shown in red. 



Energy (keV) 


Fig. 11.— Measured PSF fitting parameters as a function of 
energy. 


tive PSFs) and backgrounds, we calculated y 2 for vari¬ 
ous values of fitting parameters for each observation, and 
found a set of parameters that minimize the combined 
y 2 . An example the observed radial profile for a CygX-1 
observation (obs. ID: 10002003001) and the best-fit PSF 
model is shown in Figure [To) 

From the analysis, we found that the PSF sharpens 
- the half power diameter decreases - with energy. As 
a function of increasing energy we find a smaller core¬ 
convolution width and narrower wing, as well as a rela¬ 
tively larger amplitude for the core component (see Fig¬ 


ure 111. The sharpening saturates at ~10keV. The dif¬ 
ferences in the enclosed count fractions at different ener¬ 
gies are modest, with a maximum of ~3% between the 
lowest and the highest energy bands for a ~1' extrac¬ 
tion radius region as shown in Figure 12 However, since 
this is significant enough for small extraction regions, we 
include the energy dependent evolution into the ARF 
generation, linearly interpolating between energies. The 
additional broadening of the wings of the PSF at ener¬ 
gies below lOkeV is likely explained by scattering from 
surface contaminants. 
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TABLE 7 

Gain-Corrections (2012-vs-2015) 


Detector Number 

A Slope 

A Offset (eV) 

B Slope 

B Offset (eV) 

0 

0.9944 

-60 

0.9964 

-59 

1 

0.9928 

-66 

0.9939 

-56 

2 

0.9999 

-34 

0.9937 

-34 

3 

0.9977 

-55 

0.9963 

-7 


TABLE 8 

Cross-calibration Campaign Observations 


Observatory 

OBSID 

Exposure 

(ksec) 

OBSID 

Exposure 

(ksec) 


PKS2155-304 

3C273 

NuSTAR 

60002022002 

45 

10002020001 

244 

Chandra 

15475 

30 

14455 

30 

X MM-Newton 

0411782101 

76 

0414191001 

38.9 

Swift 

00030795108 

17.7 

00050900019 

13 

Swift 

- 

- 

00050900020 

6.9 

Suzaku 

108010010 

53.3 

107013010 

39.8 



Radius (") 


Fig. 12.— Enclosed counts fraction in different energy bands 
(top) and the differences from that in the 3—4.5 keV band (bottom). 


3.1.1. Off-Axis PSF 

The FWHM of the NuSTAR PSF does not change with 
off-axis angle, but the shape of the PSF gradually dis¬ 
torts azimuthally, and the geometrical shadowing of the 
shells causes the PSF to appear elongated, as shown in 
Figure 13 Probing the PSF at high off-axis angles is 
complicated by the lower count-rate, and necessitates a 
2-D approach where we match the contours of the model 
PSF to the actual PSF. 

We followed the same procedure (e.g., convolution, 
aspect-weighting, etc) as we did for the on-axis studies 
to produce the model PSF for the off-axis observations, 
6 > 3' (listed in Table [f]). The ray-trace PSF includes 
the effect of geometrical shadowing (vignetting of pho¬ 
tons), and successfully reproduces the observed PSF at 
the core. However, the wing component, which domi¬ 
nates at R>100", did not match with the observed con¬ 
tours at large radii. We parametrized the distortion of 
the PSF wing with an ellipticity factor and a shadow 
function, where the form of the shadow function was de- 
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Fig. 13.— Upper: Source image. Lower: Count contours (solid) 
and PSF contours (dashed) for a ~~3 / off-axis observation of the 
Cyg X-l in a 2' X 2' field. 


































11 



Fig. 14.— Observed 3-4.5 keV PSFs (top) for two Vela X-l 
separated by ~300 days and the difference with lcr statistical errors 
(bottom). The PSF at T=0 is shown in black, and PSFs measured 
at T=300 days are shown in blue and red. 


rived using the geometry of the optics: 

-F s hadow(£, y,0) = 1 - A(9)cos 2 (ta,n~ 1 (B(d)y/x)), 

F vfing (x,y,0) = e- a V* 2 + B «»y\ 

where A{9) is the shadow factor, B{9) is the ellipticity 
factor, and 6 is the off-axis angle in arcminutes. The 
shadow function was multiplied onto the wing compo¬ 
nent, which was added to the ray-trace PSF. The pa¬ 
rameters were determined using the off-axis point source 
observations to match the 2-D contours at large radii, 
and we found that A(0) = 0.0250 and B(9) = 1 — 0.0250 
made the contours match at large radii (see Figure [l3|). 

With these correction factors we generated new energy- 
dependent 2-D PSF hies for the NuSTAR CALDB at 
3-4.5, 4.5-6, 6-8, 8-12, 12-20, and 20-78 keV. 


3.1.2. PSF Stability over Time 

Changes in PSF over time can occur from out-gassing 
of the epoxy used to bond glass substrates, and/or by 
temperature gradients applied to the optics by the Sun. 
Both have been intensively studied with ground experi¬ 
ments, which showed that the PSF should be stable over 
~10 years. Out-gassing is expected to be strongest in the 
first years after l aunch , and we compared observations of 
Vela X-l ( Fiirst et al.||2014 ; Table Tsl) taken ~ 300 days 
apart to look for broadening in the PSF that might sug¬ 
gest contamination or out-gassing had occurred. 

We constructed radial prof iles f or the three observa¬ 
tions of Vela X-l (Fiirst et al.|20141 for the energy bands 
used above, and compared the radial profiles directly for 
each radial bin. We subtracted the late-time profiles from 
the reference profile (T=0, obs. ID 10002007001) in or¬ 
der to see if there was a significant change in the radial 
profiles. The radial profiles and t he d ifferences in the 
3-4.5 keV band are shown in Figure pT} We find there is 
no significant change in the radial profiles. 


4. DETECTOR GAIN CALIBRATION 


Each NuSTAR focal plane contains a calibration source 
primarily composed of 155 Eu which has several promi¬ 
nent lines at high (86.54 and 105.4 keV) and low energies 
(6.06 and 6.71 keV). The source is located on a movable 
arm that can deploy the calibration source into the field 
of view so that we can confirm the gain of the instrument 
in flight. 

We obtained two epochs of calibration source data: 
one just after observatory commissioning in June 2012 
and one after the prime mission completed in January 
2015. The initial values of the CALDB files were deter¬ 
mined using extensive pixel-by-pixel-by-grade calibration 
of the detectors on the ground, where grade r efers t o the 


charg e sharing pattern among pixels (see Kitaguchi et al. 


2011 for details). We used the first epocn of calibration 


data to fine-tune the CLC, “Charge Loss Coefficient”, and 
FPM_Gain CALDB Hies based on the in-orbit conditions 
of the instrument. 

As we have limited on-orbit data, we determine the 
CLC corrections for each detector (e.g. integrated over 
all pixels) separately for grades 0, 1, 2, 3, and 4 (corre¬ 
sponding to one- and two-pixel events) and as a block 
for grades 5 through 8 (three-pixel events) and a single 
correction for grades 9 through 20. We used the second 
epoch to introduce a time-dependent gain correction in 
the FPM_gain CALDB file that takes out an apparent 
«0.2% per year shift in the gain as well as a residual 
small shift in the zero-point offset of the energy scale 
(hereafter “slope” and “offset”, respectively) using the 
following equation: 


^new = ^° kL - Offset. (5) 

Slope 

We compare the first epoch data to a reference model 
spectrum to determine the absolute energy scale cor¬ 
rections required (if any) using Equation p] We con¬ 
structed a GEANT4 model of the calibration source dur¬ 
ing ground calibration and produced a simulated spec¬ 
trum based on the blend of Eu isotopes contained in the 
calibration source. We multiply this input spectral model 
by an absorption curve based on the attenuation ma¬ 
terial between the calibration source and the detectors 
(primarily the Beryllium entrance windows to the focal 
plane) and then convolve the result with the detector re¬ 
sponse matrix (RMF) stored in the CALDB. This results 
in a “counts” spectrum that can be directly compared to 
the observed spectrum to test for changes in the energy 
scaling. 

We fine-tune the CLC parameters by performing a fit 
over a change in the slope and the offset and minimiz¬ 
ing the resulting chi-square value. The resulting model 
fits result in large reduced chi-square values (> 5) owing 
to the large offset between the model and the observed 
spectrum at low energies (the result of an unmodeled 
spectral component likely due to the Compton scattering 
of gamma-ray photons in the detectors) and small differ¬ 
ences in the shape of the 86.54 keV line (only statistically 
significant due to the large number of source counts near 
86 keV). However, the line centroids that we recover are 
accurate across the NuSTAR science bandpass from 3 to 
79 keV (Figure [l5|. 

Iterating this analysis (e.g. applying the correction and 
then re-fitting the “new” observed spectrum to the model 
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Fig. 15.— Representative spectrum of the calibration source for FPMB Det 0 (black) compared with the spectral model of the calibration 
source convolved through the instrument response (red, see text). Left: the lowest energy lines along with the locations of the 6.07 and 
6.71 keV X-ray lines in the calibration source (vertical green lines.); Middle: The 35 to 55 keV bandpass; Right: the 55 to 88 keV bandpass 
containing the main 86.54 keV line. 





Fig. 16.— Representative spectrum of the calibration source for 2012 epoch FPMB from detector 0 (black) compared with 2015 epoch 
data from the same detector. Left: A zoomed comparison of the region near the 86.54 keV line show the 2012 data (black) and the 2015 
data before (green) and after (magenta) the time-dependent correction has been applied. Middle: A comparison of the 2012 (black) and 
the corrected 2015 data (magenta) at low-energies; Right: A comparison of the 2012 (black) and corrected 2015 (magenta) data across the 
55 to 88 keV bandpass. 


spectrum) results in residual fits with slope values typi¬ 
cally within 2 x 1CU 4 of unity and offset values within 40 
eV of zero. We note that the NuSTAR data are binned 
to a native resolution of 40 eV, so the fact that we find 
residual offset errors on that scale is likely evidence for 
some aliasing. We therefore recommend adopting sys¬ 
tematic errors of 40 eV on the offset and 2 x 10 -4 on the 
slope. This implies a systematic uncertainty of 40 eV at 
energies near Fe emission features in the spectrum and 
a systematic uncertainty of ~ 60 eV near the 67.86 keV 
line used for the analysis of 44 Ti (40 eV uncertainty in 
the offset combined with ~ 20 eV due to the uncertainty 
in the slope at 68 keV). 

Using the second epoch of calibration source data we 
can determine if there is a long-term trend in the gain 
of the instrument. We do this by comparing the 2015 
epoch data to the 2012 and adjusting the energy channel 
(PI) value and minimizing the chi-square value between 
the two distributions. Again, subtle changes in the line 
shape with time near the 86 keV line and the large num¬ 
ber of counts in that line result in a reduced chi-square 
value that is formally unacceptable (and therefore make 
it impossible to quote formal statistical errors on the fit 
parameters). However, we do find by visual inspection 
that the adjusted 2015 spectra do match the 2012 at all 
energies (Figure |l6|)). 

The corrections required to scale the 2015 up to match 
the 2012 epoch observations are found in Table [7] We 
note that, as above, once these corrections are applied 
we find residual fits with slopes within 2 x 1CU 4 of unity 


and offsets within 40 eV of zero, so we consider these the 
residual uncertainties on this fitting method. 

The time-dependent gains are implemented in the stan¬ 
dard NuSTARDAS pipeline. As of this writing we as¬ 
sume that the time-dependence in linear and can be in¬ 
terpolated between the 2012 and 2015 data sets and can 
be extrapolated beyond 2015. 

5. ABSOLUTE NORMALIZATION AND 
CROSS-CALIBRATION 

A set of observations early in the mission, which were 
individually simultaneous with Swift, Suzaku, XMM and 
Chandra, guided our decision to adjust the final absolute 
normalization by 15%, thereby setting the Crab spec¬ 
trum normalization as measured on average by NuSTAR 
to N=8.5 keV -1 cm -2 s _1 at 1 keV. This is a flat ad¬ 
justment at all energies and applied directly into the 
on-axis effective area files to affect all off-axis angles 
through the multiplication of the vignetting function. 
NuSTAR. participated in the cross-calibration campaign 
between Chandra, Suzaku, Swift and XMM-Newton per¬ 
formed on the quasar 3C273 on UT 2012 July 17 and 
the blazar PKS2155-304 on UT 2013 April 23, and we 
present here using these two sources the results of the 
cross-calibration after the normalization correction, ad¬ 
dressing the two important questions: (i) How different 
is NuSTAR with respect to the above mentioned obser¬ 
vatories? (ii) Accepting the differences, what are the 
cross-normalization constants that yield the best fit of a 
model between NuSTAR and the other observatories? 
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Fig. 17.— Contours of 1, 2 and 3 cr of the individual fits of each instrument. Top row: PKS 2155-304. Bottom rows: 3C273. 


Cross-calibration is an ongoing effort coordinated by 
the International Astronomical Consortium for High En¬ 
ergy Calibration, IACHEC^] and this analysis is a subset 
of more detailed forthcoming paper on this campaign. 
Future cross-calibration campaigns of the IACHEC will 
expand on and modify these results as the response files 
of the involved instruments may change. 

5.1. Data Reduction 

The cross-calibration observations are listed in Ta¬ 
ble [8] and were processed using the respective standard 
pipelines. 

For Chandra we used CIAO 4.6.1 and CALDB 4.6.1.1. 
For 3C 273 the data were taken with gratings configu¬ 
ration ACIS+HETG and reprocessed using the CIAO 
chandra_repro reprocessing script. We combined orders 
1-3 for the HEG and MEG arm separately, and binned 
the data at 30 counts. For PKS 2155-304 the data were 
taken with gratings configuration ACIS+LETG. We fit 
orders ml and pi simultaneously. 

For NuSTAR we used HEAsoft 6.15.1 and CALDB 
20131223. The data were processed with all standard 
settings and source counts extracted from a 30" radius 
circular region for both 3C 273 and PKS 2155-304. Back¬ 
ground was taken from the same detector. We did not 

2 http://web.mit.edu/iachec/ 


combine modules. 

For Suzaku we used CALDB: HXD (20110913), XIS 
(20140203). For 3C 273 the observation was taken in 1/4 
window mode, and we used 100 arcsec radius circular re¬ 
gions for the front-illuminates, FI, detectors (XIS0,3) and 
a 140 arcsec region for the back-illuminated, BI, detec¬ 
tor (XIS1), such that the regions were centered on the 
source, but were restricted to the operational portions 
of the detectors. For PKS 2155-304 the observation was 
also taken in 1/4 window mode. We extracted from 115" 
radius circular region for the FI and 135" from the BI. 

For Swift we used HEAsoft 6.15.1 and XRT CALDB 
2014-02-04. The data for 3C 273 were taken in ‘PHO¬ 
TON’ mode and PKS 2155-304 in ‘WINDOW TIMING’ 
mode. Both were reduced using xrtpipeline. Spectra 
for 3C 273 were extracted from an annulus region, inner 
radius 5" and outer radius 30" to correct for pileup. The 
two observations for 3C 273 were combined. 

For XMM-Newton we used SAS v. 13.5.0 with CALDB 
2014-01-31. The data for 3C273 were taken in ‘Small 
Window’ mode, and to corrected for pileup in the MOS 
we extracted counts form an annulus region with inner 
radius of 15" and outer radius of 45" . For the PN we 
extracted from a circular region of radius 45" . For 
PKS 2155-304 the observation was also taken in ‘Small 
Window’ mode, and for all instruments we extracted 
from an annular region of inner radius 10" and outer 
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radius 36". 

Due to the relative beating of the South Atlantic 
Anomaly (SAA) passages and occupation periods be¬ 
tween the low Earth orbit observatories ( NuSTAR, 
Suzaku, and Swift), we decided to forego strict simul¬ 
taneity between the observatories. Instead we compare 
NuSTAR exclusively to each instrument separately and 
truncate the observations to have matching START and 
STOP times, ignoring SAA passages and occupation pe¬ 
riods. 

For the data anal ysis we u se the XSPEC version 12.8.2 
analysis sof tware ( Arnaud 1996), fitting with cstat 
statistics (Cash 1970) and present lcr errors unless oth¬ 
erwise stated. 

5.2. Cross-calibration results 

We fitted the spectra independently in XSPEC with an 
absorbed pegged power-law (tbabsxpegpwrlw) since it 
uses as normalization the flux between two energies and 
breaks the coupling of the regular model between power- 
law index and normalization at 1 keV. In this manner 
we can distinguish between slope differences and over¬ 
all normalization offsets. The hydrogen column was held 
fixed in all cases at 1.79 xlO 20 em~ 2 for 3C273 and 1.42 
xlO 20 


.. — 2 


cm ~ for PKS2155-304 (Dickey & Lockman 1990). 


The flux normalization of the pegpwrlw was calculated 
between 3-7 keV for the soft instruments and 20-40 keV 
for the hard instruments. We limited the NuSTAR fitting 
range to match that of the comparison instrument and 
set the lower energy for the soft instruments to 3 keV, ex¬ 
cept for Chandra and Swift where it was necessary to go 


made simultaneous fits between each instrument 
respectively and FPMA and FPMB assuming 
constantxtbabsxpegpwrlw. We tied all parame¬ 
ters together and only allowed the constant between 
them to float to evaluate a potential skewing the slope 
differences might cause on the cross-calibration constant 
compared to the ratio of fluxes. The results are sum¬ 
marized in Table [lOj In all cases the constants are very 
similar to the FPMA normalized flux from Table 1 We 
do, however, strongly encourage users to first evaluate 
the parameters of individual fits and calculate fluxes 
before before applying a cross-normalization constant 
and tying fit parameters together. Because of the pos¬ 
sible variations in flux between the NuSTAR FPMs, we 
advise applying a floating cross-normalization constant, 
but with the caveat that if it takes on extreme values 
the data reductions might have to be re-examined, and 
if significant residuals occur after fitting it may be an 
indication that the source signal-to-noise is high enough 
to make inter-instrumental slope differences significant. 
Large discrepancies between FPMA and FPMB may be 
due to the interference of a gap, which will not affect the 
shape of the spectrum, but cause one module to have a 
lower flux. 

Different cross-normalization constants may occur if 
fluxes for energy ranges other than those presented here 
are used. This is due to slope differences between ob¬ 
servatories, and if simultaneity between observations is 
not strictly enforced this may also have an effect. In 
general, however, typical values found for other sources 
than those discussed above are within 10% as demon¬ 


strated in Walton et al. (2014 1 for Holmberg IX X-l 


may affect the slope measurement, it does not affect the 

(XMM-Newton and Suzaku F, Brenneman et al. 

(20141 

flux calculation. Lhe upper range was 8keV tor Chan- 

lor IC 4329A ( Suzaku) and Balokovic et al. (20 

M for 

dra, 8keV for Swift, 9keV for XMM-Newton and 9kcV 

NGC 424, NGC 1320, and IG 2550 (. XMM-Newton and 


for Suzaku. 

The resul ts a re summarized in Table[9]for both sources, 
and Figure [TT! shows the 1, 2 and 3a contours of the in¬ 
dividual fits of the instruments. In all cases there is con¬ 
sistency between the measured slopes when considering 
the 3a confidence contours, but obvious slope differences 
do exist, most notably between NuSTAR and XMM- 
newton/ MOS instruments and Swift/XRT in 3C273. It 
is possible that this discrepancy is in part due to the 
brightness of 3C 273, which caused pileup in both instru¬ 
ments. While excision of the piled-up region is common 
practice, the errors in the PSF wings of the instruments 
start to play a role and may have skewed the results. 

The ratio of fluxes are shown relative to FPMA and 
NuSTAR was intentionally calibrated to be roughly in 
the middle of the spread between instruments with Chan¬ 
dra gratings and Suzaku /HXD yielding the highest fluxes 
(Rfpma ~ 1.10) and Swift /XRT generally the lowest 
flux es b etween 3-7keV (Rfpma ~ 0.95). As discussed 
in § |2.4| we expect to see differences between FPMA and 
FPMB due to the restrictions on the absolute optical axis 
knowledge and while for PKS2155-304 the two modules 
yielded the same flux, during 3C 273 there was an offset of 
3%. It is well known that in gener a l XMM -Newton/ pn is 
low compared to MOS ( Read et al.|2014| ), but for 3C 273 


Swift). 

We have in the above not included the ^3% uncer¬ 
tainty on the NuSTAR fluxes, which are due to the con¬ 
straints on absolute knowle dge o f the optical axis loca¬ 
tion as discussed in Section |2~4l Finally we stress that 
since both targets are variable and observations were of 
different lengths between instruments, the non-NuSTAR 
instrument fluxes should not be compared against each 
oth er, only w ith re spect to NuSTAR. A forthcoming pa¬ 
per Madsen (prep) of the IACHEC consortium on these 


same cross-calibration campaigns, will discuss the cross¬ 
calibration between all instruments in far greater detail. 

6. TIMING CALIBRATION 
6.1. Estimated Performance 
A basic de scription of the NuST AR timi ng system was 


the spectrum appears to be unusually low for reasons 
that are not understood. This and the slope differences 


will be investigated in Madsen (prep). 
To investigate cross-norma' 


ization constants, we 


provid ed in Harrison et al. (2013) and Baclietti et al. 
(2015). We elaborate upon the system here and describe 
efforts to characterize the timing performance. The per¬ 
formance of the components is summarized in Table |11[ 
and described in more detail below. Known biases in the 
system are shown in one column, and are typically re¬ 
moved by the standard analysis. The unmodeled terms 
are shown in the second column, and are either unknown 
at the time of writing, or are stochastic. Based on this 
analysis, we expect the corrected event time stamps to 
be accurate at the 3 ms level (95%). 

The NuSTAR on-board time reference is a 
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TABLE 9 

Cross-calibration tbabs x pegpwrlw 


TABLE 10 

Cross-calibration constant x tbabs x pegpwrlw 


Instrument 

10 

Flux 3-7 keV 
(10 -12 erg cm 2 s 1 ) 

Rfpma[j 

PKS 2155-304 

FPMA 

2.81 ± 0.13 

5.9 ± 0.2 

1.0 

FPMB 

2.77 ± 0.13 

5.8 ± 0.2 

0.98 ±0.05 

ACIS LETGS 

2.86 ±0.10 

6.6 ± 0.3 

1.12 ±0.06 

FPMA 

2.78 ± 0.07 

5.3 ± 0.1 

1.0 

FPMB 

2.75 ± 0.07 

5.3 ± 0.1 

1.0 ± 0.03 

XIS0 

2.8 ± 0.09 

5.2 ±0.1 

0.98 ± 0.03 

XIS1 

2.84 ± 0.09 

5.1 ±0.1 

0.96 ± 0.03 

XIS3 

2.87 ± 0.08 

5.2 ±0.1 

0.98 ± 0.03 

FPMA 

2.82 ± 0.1 

5.5 ± 0.1 

1.0 

FPMB 

2.74 ± 0.1 

5.3 ± 0.1 

0.96 ± 0.03 

XRT 

2.88 ± 0.1 

5.2 ± 0.3 

0.95 ± 0.06 

FPMA 

2.8 ± 0.08 

5.6 ± 0.1 

1.0 

FPMB 

2.8 ± 0.08 

5.6 ± 0.1 

1.0 ± 0.03 

MOS1 

2.71 ±0.04 

5.8 ±0.1 

1.04 ± 0.03 

MOS2 

2.76 ±0.04 

5.6 ±0.1 

1.0 ± 0.03 

pn 

2.78 ±0.04 

5.5 ±0.1 

0.98 ± 0.03 


3C273 


FPMA 

1.52 

± 

0.09 

40.3 

± 

0.6 

1.0 

FPMB 

1.68 

± 

0.09 

41.3 

± 

0.7 

1.02 ±0.02 

ACIS HEGTS 

1.52 

±0.05 

44.6 

± 

1.0 

1.10 ± 0.03 

FPMA 

1.66 

± 

0.02 

40.9 

± 

0.3 

1.0 

FPMB 

1.66 

± 

0.02 

42.5 

± 

0.3 

1.04 ± 0.01 

XIS0 

1.64 

± 

0.03 

40.6 

± 

0.3 

0.99 ± 0.01 

XIS1 

1.69 

± 

0.03 

40.4 

± 

0.3 

0.99 ± 0.01 

XIS3 

1.66 

± 

0.03 

41.6 

± 

0.3 

1.02 ± 0.01 

FPMA 

1.65 

± 

0.03 

40.9 

± 

0.4 

1.0 

FPMB 

1.66 

± 

0.03 

42.4 

± 

0.4 

1.04 ± 0.01 

XRT 

1.55 

± 

0.06 

42.4 

± 

1.2 

1.04 ± 0.03 

FPMA 

1.63 

± 

0.04 

39.8 

± 

0.6 

1.0 

FPMB 

1.68 

± 

0.04 

41.3 

± 

0.6 

1.04 ± 0.02 

MOS1 

1.47 

±0.04 

40.7 

±0.4 

1.02 ± 0.02 

MOS2 

1.46 

±0.04 

39.0 

±0.4 

0.98 ± 0.02 

pn 

1.59 

±0.02 

37.0 

±0.2 

0.93 ± 0.01 


Instrument 

r 

Flux 20-40 keV 
(10 -12 erg cm 2 s 1 ) 

Relative 

Flux 

FPMA 

1.73 ± 0.06 

60.4 ± 1.4 

1.0 

FPMB 

1.83 ± 0.06 

60.6 ± 1.4 

1.0 ± 0.03 

HXD 

1.75 ± 0.02 

68 ±3 

1.12 ± 0.03 


individual fits to each instrument pair. 

b Ratio of flux relative to FPMA, Flux y/F1uxfpma 


Instrument 

rQ 

Flux 3-7 keV 
(10 -12 erg cm 2 s 1 ) 

Constanipl 

PKS 2155-304 


ACIS LETGS 

2.83 ±0.06 

5.9 ± 0.2 

1.15 ± 0.06 

XIS0 

2.80 ± 0.05 

5.5 ±0.1 

0.95 ±0.03 

XIS1 

2.81 ± 0.05 

5.5 ±0.1 

0.93 ±0.03 

XIS3 

2.82 ± 0.04 

5.5 ±0.1 

0.94 ±0.03 

XRT 

2.83 ± 0.05 

5.4 ± 0.2 

1.00 ±0.05 

MOS1 

2.76 ±0.02 

5.6 ±0.1 

1.03 ±0.03 

MOS2 

2.77 ±0.02 

5.6 ±0.1 

1.00 ±0.03 

pn 

2.74 ±0.02 

5.6 ±0.1 

0.98 ±0.03 


3C273 


ACIS HETGS 

1.57 ±0.05 

40.0 ± 0.6 

1.10 ± 0.04 

XIS0 

1.65 ± 0.01 

40.9 ±0.3 

0.99 ±0.01 

XIS1 

1.67 ± 0.01 

41.1 ±0.3 

0.99 ±0.01 

XIS3 

1.66 ± 0.01 

41.0 ±0.3 

0.01 ±0.01 

XRT 

1.63 ± 0.03 

40.3 ± 0.5 

0.99 ±0.03 

MOS1 

1.63 ±0.03 

39 ±0.5 

1.05 ±0.02 

MOS2 

1.62 ±0.02 

39.7 ±0.5 

1.06 ±0.02 

pn 

1.65 ±0.02 

40.0 ±0.5 

0.93 ±0.01 


Instrument 

r 

Flux 20-40 keV 
(10 -12 erg cm 2 s 1 ) 

Constant 

HXD 

1.77 ± 0.04 

60 ±1 

1.13 ±0.04 


Simultaneous fits of each instrument with NuSTAR, only Con¬ 
stant allowed to float 

b ConstantppMA — 1. 


which allows measurement of the spacecraft clock offset 
and trends. The offset is measured essentially by per¬ 
forming one-way ranging of telemetry frames set 1 sec¬ 
ond intervals. The time of transmission as recorded on 
the spacecraft and the time of reception on the ground as 
recorded by the ground station are used to estimate the 
clock offset. During downlink contacts, the rising edge 
of the PPS signal initiates the downlink of a telemetry 
frame. There is a delay of at least 352.75 ms before this 
packet is transmitted to the ground. During times when 
the downlink system is under heavy load, this delay can 
be larger, so clock offset data are filtered to consider only 
idle downlink times. One-way range delays are modeled 
using the known spacecraft ephemeris, which is typically 
accurate to ~30 km (^100 lt-/zs). The ground stations 


temperature-compensated crystal oscillator with a 
nominal frequency of 24 MHz. The frequency stability 
of this oscillator is specified to be ±4 ppm over the 
full operating temperature range. On-board time 
synchronization of subsystems is done by broadcasting 
a one pulse-per-second (PPS) electronic signal. The 
propagation delay of this signal is expected to be less 
than a few microseconds. The instrument focal plane 
modules maintain their own internal high resolution 
timers: seconds are marked using the PPS signal, and 
sub-seconds are marked internally with 2 /is precision. 
The delay between when a photon interacts with the 
detector and when it is time-stamped by the electronics 
is expected to be a few microseconds, but is certainly 
less than 0.1msec. After the spacecraft clock timestamp 
is applied to an event, it is telemetered to the ground. 
Science event data is archived with this timestamp. 

The spacecraft clock is routinely correlated with 
ground time references during most ground contacts, 


maintain their own time standards referenced to UTC via 
GPS, which are typically accurate at the few nanosecond 
level. However, USN Hawaii and KSAT Singapore sta¬ 
tions only record receive times with 1 ms precision. The 
RF propagation and Malindi ground station portions of 
the system have been validated by analysis of Swift data 
to better than 0.1ms ( Cusumano et al.]|2012 ). 

The clock offset trend is examined on the ground. The 
NuSTAR system is required to maintain on-board time 
to within 100 ms of UTC, but typically the offset is main¬ 
tained to less than 30 ms. As the spacecraft clock’s envi- 
roment changes, the true frequency also changes. Mission 
operations staff can adjust the clock divisor setting from 
the nominal 24 MHz value in order to steer the clock 
offset towards zero. The clock divisor is the “effective” 
frequency, or oscillator ticks required to produce a PPS 
pulse. The mission-average clock divisor offset is ~14 
ppm, with a variation of about ±0.2 ppm (95%). Clock 
divisor adjustments are visible in both the science event 
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TABLE 11 

NuSTAR Expected Timing Performance 


Error Term 

Modeled Bias 
(ms) 

Unmodeled 

(ms) 

Spacecraft 1 PPS time transfer 

— 

<0.1 

Instrument time stamp 

— 

<0.1 

Polynomial model 

<100 

<2 

Spacecraft 1 PPS downlink delay 

352.75 

<!0 

Range delay 

<10 

<0.1 

Ground station time stamp 

— 

<°. .ip 


a During idle downlink periods 
b Malindi ground station 
C USN Hawaii ground station 
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Fig. 18. — Comparison of the TOAs of the Crab and the ex¬ 
pected arrival times calculated through the Jodrell Bank Monthly 
Ephemeris in all available NuSTAR observations. The historical 
residuals are about ±3 ms. The inset shows a detail of the October 
2013 observation. This plot shows a clear shift by about 0.4±0.1 ms 
of the measured TOAs, originated by clock frequency variations, 
on orbital timescales. Dashed lines indicate the orbital timescale 
for visual purposes. 


timestamps and the clock offset measurements, and are 
included in the clock modeling as described below. 

We model the clock offset trend as a piece-wise contin¬ 
uous polynomial function. Control points are inserted at 
points in time when the clock divisor was adjusted, and 
also where residuals demand extra compliance. The coef¬ 
ficients of the function are recorded in a clock correction 
file, which is available to the NuSTAR community via the 
CALDB system. After modeling, the residuals are less 
than 2 ms (95%), although in cases where the clock mea¬ 
surements are sparse, the retrieval error may be larger 
than this. This file can be used by the general-purpose 
FTOOL barycorr in order to correct event times to the 
solar system barycenter time system, or TDB. As the 
NuSTAR orbit is known to approximately 30 km, light 
time errors would be less than 0.1 ms. 

6.2. Measured Performance 

We used observations of known X-ray pulsars in or¬ 
der to characterize the actual performance of the sys¬ 
tem. This includes absolute timing, as well as long term 
stability of the timing measurement system. We ana¬ 
lyzed the extensive observations of the Crab pulsar de¬ 
scribed above. Observation of a single pulsar allows one 


to track the mission-long stability of the clock, mod¬ 
ulo the ~33ms pulse period. We also observed PSR 
B1509—58, with a longer period of 151 ms to resolve some 
of this ambiguity. We describe the results of these obser¬ 
vations below. We refer observations of known pulsars 
to their absolute radio ephemerides. 

We used all available NuSTAR observations of the 
Crab pulsar (Table [I]), and extracted pulse profi les i n 
the 3-79keV range. We used the PRESTC^] (Ransom 


200l|) suite of programs to obtain an average pulsed pro- 
file and to calculate the time-of-arrival (TOA) of the 
pulses on time scales of ~ 5 minutes. We then com¬ 
pared these TOAs to the relevant monthly timing so¬ 
lution for th e Crab produced by the Jodrell Bank radio 
telescope F|(Lyne et ah 1993) using the Tempo2 software^] 
(Edwards et al. 2(J(J6|). The radio ephemeris is referrea 
to well-known terrestrial time standards. There is an 
offset between radio and X-ray pulse profiles of ~300 /is, 
which we do not consider further in this analysis because 
it is negligible in light of the current NuSTAR timing 
performance. 

The results are shown in Figure [18] Over the whole 
history of NuSTAR Crab observations, the residuals of 
these TOAs show a trend between —2 ms and +3 ms. 
This level of residual is indeed mostly inside the ± ~ 3 ms 
uncertainty range estimated above. 

In addition to the long term residual trend, sh orte r 
term trends are also visible (see inset of Figure 18). 
These trends repeat on the time scale of the spacecraft 
orbital period of about 5400s (0.07 day), and have an 
amplitude of about 400 /its. This trend is clearly related 
to the temperature history of the spacecraft clock oscil¬ 
lator. Oscillator instability on these time scales hinders 
the detection of pulsars with period shorter than 3 ms, 
because the pulse peak will be smeared over a spacecraft 
orbital period. 

We also measured the TOAs of another pulsar, PSR 
B1509—58, with a pulse period of ~ 151ms and verified 
that the TOAs were aligned with those measured by Swift 
in SW timing mode to 1 ± 2 ms (details can be found in 


Mori et al. (2014). 

In summary, using two pulsars with known absolute 
ephemerides, the NuSTAR timing system is performing 
at the accuracy level of 3 ms over the entire mission to 
date, which is about the expected level of timing uncer¬ 
tainty based on design (Table [TTj) . 

Present research indicates that the spacecraft clock fre¬ 
quency variations are driven almost entirely by its ther¬ 
mal environment. As accurate temperature information 
about the environment is known, it may be possible to 
make a more accurate clock variation model. This will 
be the subject of future work. 


7. SUMMARY 

The NuSTAR observatory is well understood, cali¬ 
brated, and operating within requirements. We have 
calibrated the instrument effective area responses against 
the phase-averaged Crab nebula + pulsar out to T as¬ 
suming a canonical spectrum of P = 2.1 and normal¬ 
ization N = 8.5. The scatter in fit values of the Crab 


3 http://www.cv.nrao.edu/ sransom/presto/ 

4 http://www.jb.man.ac.uk/pulsar/crab.html 

5 http://www.atnf.csiro.au/research/pulsar/tempo2/ 
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data-set after adjustment are summarized in Tableland 
are a measurement of the repeatability between consec¬ 
utive observation of the same source with varying off- 
axis angle (< 3') and position on the detector. For the 
Crab ensemble we find T = 2.1 ± 0.01 and normalization 
N = 8.5 ± 0.3. The flux differences between FPMA and 
FPMB are of the order 0-5% and understood to be due 
to mast motions and the uncertainty in the true optical 
axis location, which is known to ~ 30". The detectors 
have between 2012 and 2015 shown a gain shift in slope 
and offset (summarized in Tableftl and the estimated un¬ 
certainty on the energy is 60eVTCross-calibration cam¬ 
paigns with Chandra , Swift, Suzaku and XMM-Newton 
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